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Introduction. The article considers the dynamics of air pollution by stationary sources in the Russian Federation from 
1998 to 2016. Harmful emissions into the atmosphere cause great harm to all living organisms. As a result, the life 
expectancy of the population is significantly reduced. Therefore, the assessment of pollution volumes and subsequent 
measures for the protection of atmospheric air are priority tasks of our time. 

Problem statement. The objectives of the study are to analyze the dynamics of pollution, build a mathematical model of 
this process, and implement a forecast for a five-year period. 

Theoretical part. Data for the work is taken from the official statistical book. Microsoft Excel and StatSoftStatistica 
computer technologies are used for calculations. 

Conclusion. Based on the analysis, an adequate mathematical model is constructed, which may be of interest for 
predicting the anthropogenic impact on the environment. 
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Introduction. Air pollution causes great harm not only to living organisms, but also to the hydrosphere, soil 
and vegetation cover, buildings, structures and other objects. People who live in places with polluted air often have 
diseases such as allergies, cancer, and various lung diseases. As a result, life expectancy is significantly reduced. 
Therefore, the protection of atmospheric air is one of the priorities of our time. 

Main sources of atmospheric pollution in Russia are: industry, transport, utilities and agriculture. The level of 
air pollution depends, as a rule, on the degree of urbanization and industrial development of the territory. 

The paper considers the amount of harmful substances (hereinafter — pollution) released by stationary sources 
into the atmosphere in 1998-2016. The task of the authors is to analyze the dynamics of changes in pollution, to build a 
mathematical model of this process and to forecast it for a five-year period. The source data is taken from the book 
"Russia in numbers" [1]. 

Methodology. The values of atmospheric air pollution are presented in Table 1, and Fig. 1 shows their 
graphical representation. 



































Table 1 
Air pollution values in the Russian Federation 

Year Pollution, million tons Year Pollution, milliontons 
1998 18,7 2008 20,1 

1999 18,5 2009 19 

2000 18,8 2010 19,1 

2001 19,1 2011 19,2 

2002 19,5 2012 19,6 

2003 19,8 2013 18,4 

2004 20,5 2014 17,5 

2005 20,4 2015 17,3 

2006 20,6 2016 17,3 

2007 20,6 - - 
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Fig. 1. Dynamics of changes in air pollution in the Russian Federation 


The initial data show fluctuations in pollution against the downward trend, so the model dependence of 
pollution on time is a linear combination of monotonically decreasing and oscillatory functions. The monotonically 
decreasing function is an exponent, and the oscillating function is a sinusoid [2-4]: 


= 21 
X = A+ Be~CE-18) + D «sin [Fe + P| (1) 


whereX — the mathematical expectation of the pollution value; t — the year; 4, B, C, D, E, F — the parameters of the 
exponential-harmonic function. 

Since there are many adjustable parameters in function (1): three parameters in the exponent (constant, 
amplitude, and time) and three parameters in the sine wave (amplitude, frequency, and shift), and the initial data is 
small, the process of model identification was divided into two stages. 

At the first stage, ignoring the fluctuations, the exponential dependence was identified: 

Fest Re tee. 
The coefficients A, B and C are found by solving the optimization problem: 


xe -Y(t)) > min, (2) 


L 
wheret; — the value of the argument (time); Yi— the actual value of pollution: 


Y¥@j)=A+ Be i). 
At the second stage the oscillatory part is considered: 
= 21 
Z= D:sin [Fe +F)] 
The coefficients D, E and F are found by solving the optimization problem: 
YG) — % = VE)? + min, 3) 
i 
wheret; — the value of the argument (time): 
21 
Z(t;) =D «sin [= (t; + F)| 
Y- the actual value of pollution 
Y(t;) = A + Be~Ci-1998), 
Results.To find unknown coefficients, we used the functions of the Microsoft Excel solution search package 
[5-7]. As a result of solving the optimization problem (2) by the least squares method [8], the coefficients were found: A 
= 15,4; B = 4,378; C = 0,017. As a result of solving the optimization problem (3), the coefficients were found: D = 


1,253; E = 19,287; F = -2001,84. 
The final type of the dependency (1): 


X = 15,4 + 4,378 - e~ 9017-1998) + 1,.253- sin (t — 2001,84)}. (4) 


21 
lios87 


Table 2 shows [9] the characteristics obtained from the model (4). 
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Table 2 
Characteristicsof themodel 
: Estimated os Square of the 
Year Acta pollution, pollution, asbeolntedevial Relativedeviation, % en deviations, 
million tons faa on, milliontons 
milliontons % 
1998 18,7 18,587 0,113 0,603 0,363 
1999 18,5 18,700 0,200 1,083 1,173 
2000 18,8 18,920 0,120 0,639 0,408 
2001 19,1 19,216 0,116 0,605 0,366 
2002 19,5 19,548 0,048 0,247 0,061 
2003 19,8 19,875 0,075 0,379 0,144 
2004 20,5 20,155 0,345 1,684 2,836 
2005 20,4 20,350 0,050 0,243 0,059 
2006 20,6 20,434 0,166 0,805 0,648 
2007 20,6 20,390 0,210 1,017 1,035 
2008 20,1 20,217 0,117 0,581 0,337 
2009 19 19,925 0,925 4,866 23,678 
2010 19,1 19,538 0,438 2,293 5,258 
2011 19,2 19,091 0,109 0,567 0,322 
2012 19,6 18,624 0,976 4,977 24,774 
2013 18,4 18,181 0,219 1,191 1,419 
2014 17,5 17,801 0,301 1,718 2,951 
2015 17,3 17,518 0,218 1,259 1,585 
2016 17,3 17,356 0,056 0,323 0,105 
The following indicators are calculated: 
e the correlation coefficient is 0,936; 
e the average value of the absolute deviation is 0,253 million tons; 
e the average value of the relative deviation is 1,320 %; 
e the maximum relative deviation is 4,977 %; 
e the standard deviation is 1,885 %; 
e the model's adequacy is 1,889 %. 
Fig. 2 shows the atmospheric pollution values obtained from the model (4). 
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Fig. 2. Dynamics of changes in the atmospheric pollution in the Russian Federation: 


row | — actual data; row 2 — calculated data 


The forecast for the period 2017-2021 was also obtained using formula (4) [10]. Table 3 presents the forecast 
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Table 3 
Comparison of actual and predicted data on pollution 
Year Real value, milliontons Predictedvalue, milliontons 
2017 17,5 17,326 
2018 17,1 17,426 
2019 - 17,639 
2020 - 17,937 
2021 - 18,284 

















The forecast error is 0,994 % for 2017, and for 2018 it is 1,906 %. 

Conclusion: 

e the dynamics of air pollution from stationary sources in the Russian Federation in 1998-2016 is analyzed; 

e arelationship was found that reflects the change in the amount of harmful emissions over the specified 
period; 

e an adequate mathematical model has been constructed that may be of interest for predicting anthropogenic 
impacts on the environment and can be used as one of the tools for developing alternative models for such forecasting; 

e the forecast of the dynamics of changes in atmospheric pollution in the Russian Federation for a five-year 
period has been implemented. 
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